Understanding virus and microbial evolution in wildlife through meta-transcriptomics

Wildlife harbors a substantial and largely undocumented diversity of RNA viruses and microbial life forms. RNA viruses and microbes are also arguably the most diverse and dynamic entities on Earth. Despite their evident importance, there are major limitations in our knowledge of the diversity, ecology, and evolution of RNA viruses and microbial communities. These gaps stem from a variety of factors, including biased sampling and the difficulty in accurately identifying highly divergent sequences through sequence similarity-based analyses alone. The implementation of meta-transcriptomic sequencing has greatly contributed to narrowing this gap. In particular, the rapid increase in the number of newly described RNA viruses over the last decade provides a glimpse of the remarkable diversity within the RNA virosphere. The central goal in this thesis was to determine the diversity of RNA viruses associated with wildlife, particularly in an Australian context. To this end I exploited cutting-edge meta-transcriptomic and bioinformatic approaches to reveal the RNA virus diversity within diverse animal taxa, tissues, and environments, with a special focus on the highly divergent "dark matter" of the virome that has largely been refractory to sequence analysis. Similarly, I used these approaches to detect targeted common microbes circulating in vertebrate and invertebrate fauna. Another important goal was to assess the diversity of RNA viruses and microbes as a cornerstone within a new eco-evolutionary framework. By doing so, this thesis encompasses multiple disciplines including virus discovery, viral host-range distributions, microbial-virus and host–parasite interactions, phylogenetic analysis, and pathogen surveillance. In sum, the research presented in this thesis expands the known RNA virosphere as well as the detection and surveillance of targeted microbes in wildlife, providing new insights into the diversity, evolution, and ecology of these agents in nature

A divergent articulavirus in an Australian gecko identified using meta-transcriptomics and protein structure comparisons

The discovery of highly divergent RNA viruses is compromised by their limited sequence similarity to known viruses. Evolutionary information obtained from protein structural modelling offers a powerful approach to detect distantly related viruses based on the conservation of tertiary structures in key proteins such as the RNA-dependent RNA polymerase (RdRp). We utilised a template-based approach for protein structure prediction from amino acid sequences to identify distant evolutionary relationships among viruses detected in meta-transcriptomic sequencing data from Australian wildlife. The best predicted protein structural model was compared with the results of similarity searches against protein databases. Using this combination of meta-transcriptomics and protein structure prediction we identified the RdRp (PB1) gene segment of a divergent negative-sense RNA virus, denoted Lauta virus (LTAV), in a native Australian gecko (Gehyra lauta). The presence of this virus was confirmed by PCR and Sanger sequencing. Phylogenetic analysis revealed that Lauta virus likely represents a newly described genus within the family Amnoonviridae, order Articulavirales, that is most closely related to the fish virus Tilapia tilapinevirus (TiLV). These findings provide important insights into the evolution of negative-sense RNA viruses and structural conservation of the viral replicase among members of the order Articulavirales.This research was funded by the Australian Research Council, grant number FL170100022

Evolutionary history and spatio-temporal dynamics of dengue virus serotypes in an endemic region of Colombia.

Dengue is a prevalent disease in Colombia and all dengue virus serotypes (DENV-1 to -4) co-circulate in the country since 2001. However, the relative impact of gene flow and local diversification on epidemic dynamics is unknown due to heterogeneous sampling and lack of sufficient genetic data. The region of Santander is one of the areas with the highest incidence of dengue in Colombia. To provide a better understanding of the epidemiology of dengue, we inferred DENV population dynamics using samples collected between 1998 and 2015. We used Bayesian phylogenetic analysis and included 143 new envelope gene sequences from Colombia, mainly from the region of Santander, and 235 published sequences from representative countries in the Americas. We documented one single genotype for each serotype but multiple introductions. Whereas the majority of DENV-1, DENV-2, and DENV-4 strains fell into one single lineage, DENV-3 strains fell into two distinct lineages that co-circulated. The inferred times to the most recent common ancestors for the most recent clades of DENV-1, DENV-2, and DENV-4 fell between 1977 and 1987, and for DENV-3 was around 1995. Demographic reconstructions suggested a gradual increase in viral diversity over time. A phylogeographical analysis underscored that Colombia mainly receives viral lineages and a significant diffusion route between Colombia and Venezuela. Our findings contribute to a better understanding of the viral diversity and dengue epidemiology in Colombia

Faecal virome of the Australian grey-headed flying fox from urban/ suburban environments contains novel coronaviruses, retroviruses and sapoviruses

Bats are important reservoirs for viruses of public health and veterinary concern. Virus studies in Australian bats usually target the families Paramyxoviridae, Coronaviridae and Rhabdoviridae, with little known about their overall virome composition. We used metatranscriptomic sequencing to characterise the faecal virome of greyheaded flying foxes from three colonies in urban/suburban locations from two Australian states. We identified viruses from three mammalian-infecting (Coronaviridae, Caliciviridae, Retroviridae) and one possible mammalianinfecting (Birnaviridae) family. Of particular interest were a novel bat betacoronavirus (subgenus Nobecovirus) and a novel bat sapovirus (Caliciviridae), the first identified in Australian bats, as well as a potentially exogenous retrovirus. The novel betacoronavirus was detected in two sampling locations 1375 km apart and falls in a viral lineage likely with a long association with bats. This study highlights the utility of unbiased sequencing of faecal samples for identifying novel viruses and revealing broad-scale patterns of virus ecology and evolution

Intra-host genetic diversity of dengue virus type 4 strains from the municipality of Guarujá,

A caracterização da variabilidade genética intra-hospedeiro do vírus da dengue (DENV) é fundamental para a compreensão de sua evolução e dinâmica populacional no contexto atual como um importante patógeno viral humano. A diversidade viral acumulada em hospedeiros infectados influencia diretamente em outros aspectos como patogênese, transmissão e imunidade do hospedeiro. Contudo, apesar de existirem vários estudos sobre a diversidade genética intra-hospedeiro em dengue, nenhum tem sido relatado para DENV-4 até o presente momento. O ressurgimento e disseminação desse sorotipo foi associado com a sua co-circulação e o substituição dos sorotipos 1, 2 e 3 durante os recentes surtos no município do Guarujá-SP. Com base neste quadro epidemiológico, este estudo visou identificar a variação genética intra-hospedeiro do DENV-4 em amostras coletadas durante o surto de 2013, utilizando tecnologias de sequenciamento de nova geração. Portanto, nós caracterizamos a variabilidade genética de DENV-4 em diferentes níveis e as forcas evolutivas que afetam essa diversidade. Em adição, foram explorados os principais eventos na transmissão das variantes de DENV-4 identificadas. Nossos resultados revelaram uma baixa diversidade genética intra-hospedeiro para DENV-4. No entanto, mutação e pressões seletivas foram mecanismos importantes na variabilidade genética do vírus. A nível populacional, as variantes estão sujeitas á seleção natural negativa, não obstante identificamos seleção positiva atuando sob sítios específicos. Nenhuma evidência de recombinação foi detectada. Além disso, contra-intuitivamente, variantes de baixa frequência estão sendo transmitidas e contribuindo para diversidade genética do DENV-4 circulando em Guarujá. Nossos resultados fornecem novas evidências potencialmente úteis para futuros trabalhos focados em infecções mistas, escape imunológico, assim como o espalhamento e diversificação viral.Characterizing intra-host genetic variability in dengue virus (DENV) virus is paramount for understanding its evolution and population dynamics in the context of its current status as a major human viral pathogen. The extent to which viral diversity accrues in infected host influences aspects such as pathogenesis, transmission, and host immunity. Although there are several studies about intra-host genetic diversity in dengue, so far nothing has been revealed about DENV-4. In the Guarujá municipality in the State of São Paulo, the reemergence and spread of this serotype was associated with its co-circulation and the displacement of serotypes 1, 2 and 3 during recent outbreaks. Based on this epidemiological framework, we seek to identify the intra-host genetic variation of DENV-4 strains from samples collected during the 2013 outbreak by using deep sequencing technologies. We characterized the genetic variability of DENV-4 at different levels, and the forces shaping this diversity. Likewise, we explored major transmission events among DENV-4 variants. Our results revealed a low intra-host genetic diversity for DENV-4. However, we found selective and mutational pressures contributing to genetic diversity, while recombination did not seem play an important role. We further identified purifying selection at population level but sites subject to potential diversifying selection. Additionally, we observed low frequency haplotypes being transmitted among hosts and contributing to the viral diversity of DENV-4 circulating in Guarujá. Our findings provide preliminary insights for future studies in mixed infections, drug resistance, virus variant spread and immune scape. This study is the first effort to investigate the intra-host diversity of DENV-4

Meta-transcriptomics reveals potential virus transfer between Aedes communis mosquitoes and their parasitic water mites

Arthropods harbor a largely undocumented diversity of RNA viruses. Some arthropods, like mosquitoes, can transmit viruses to vertebrates but are themselves parasitized by other arthropod species, such as mites. Very little is known about the viruses of these ectoparasites and how they move through the host-parasite relationship. To address this, we determined the virome of both mosquitoes and the mites that feed on them. The mosquito Aedes communis is an abundant and widely distributed species in Sweden, in northern Europe. These dipterans are commonly parasitized by water mite larvae (Trombidiformes: Mideopsidae) that are hypothesized to impose negative selection pressures on the mosquito by reducing fitness. In turn, viruses are dual-host agents in the mosquito-mite interaction. We determined the RNA virus diversity of mite-free and mite-detached mosquitoes, as well as their parasitic mites, using meta-transcriptomic sequencing. Our results revealed an extensive RNA virus diversity in both mites and mosquitoes, including thirty-seven putative novel RNA viruses that cover a wide taxonomic range. Notably, a high proportion of viruses (20/37) were shared between mites and mosquitoes, while a limited number of viruses were present in a single host. Comparisons of virus composition and abundance suggest potential virus transfer between mosquitoes and mites during their symbiotic interaction. These findings shed light on virome diversity and ecology in the context of arthropod host-parasite-virus relationships

Substantial viral and bacterial diversity at the bat-tick interface

Ticks harbour a high diversity of viruses, bacteria and protozoa. The soft tick Carios vespertilionis (Argasidae) is a common ectoparasite of bats in the Palearctic region and is suspected to be vector and reservoir of viruses and other microbial species in bat populations, some of which may act as zoonotic agents for human disease. The Soprano pipistrelle (Pipistrellus pygmaeus, Vespertilionidae) is widely distributed in Europe, where it can be found inside or close to human habitation. We used metatranscriptomic sequencing to determine the RNA virome and common microbiota in blood- fed C. vespertilionis ticks collected from a Soprano pipistrelle bat roosting site in south- central Sweden. Our analyses identified 16 viruses from 11 virus families, of which 15 viruses were novel. For the first time in Sweden we identified Issuk- Kul virus, a zoonotic arthropod- borne virus previously associated with outbreaks of acute febrile illness in humans. Probable bat- associated and tick- borne viruses were classified within the families Nairoviridae, Caliciviridae and Hepeviridae, while other invertebrate- associated viruses included members of the Dicistroviridae, Iflaviridae, Nodaviridae, Partitiviridae, Permutotetraviridae, Polycipiviridae and Solemoviridae. Similarly, we found abundant bacteria in C. vespertilionis, including genera with known tick- borne bacteria, such as Coxiella spp. and Rickettsia spp. These findings demonstrate the remarkable diversity of RNA viruses and bacteria present in C. vespertilionis and highlight the importance of bat- associated ectoparasite surveillance as an effective and non- invasive means to track viruses and bacteria circulating in bats and ticks

Origin, tempo, and mode of the spread of DENV-4 Genotype IIB across the state of São Paulo, Brazil during the 2012-2013 outbreak

BACKGROUND Dengue virus type 4 (DENV-4) was first reported in Brazil in 1982 and since then no more cases were detected again in Brazil until 2010, when the virus was reintroduced. Over the following years, the virus spread to several Brazilian states and resulted in about 1,400,000 dengue cases, in 2013. The largest number of cases were documented in the Southeast macro-region. OBJECTIVES To determine the phylogeography of DENV-4 Genotype IIB strains isolated during the epidemics in 2012-2013 in São Paulo, Brazil, we aimed to contextualise the contribution of viruses sampled in different localities across the overall movement of DENV-4 in Brazil. METHODS Based on the envelope gene sequences retrieved from GenBank, we employed a Bayesian phylogeographic approach to assess the spatiotemporal dynamics of DENV-4 Genotype IIB in São Paulo, Brazil. FINDINGS The dispersal dynamics of DENV-4 Genotype IIB in Brazil indicated Rio de Janeiro and Mato Grosso states as the most likely routes toward São Paulo before the 2012-2013 outbreak. Likewise, Guarujá and São José do Rio Preto facilitated viral spread and transmission to other localities in the South and Southeast macro-regions in Brazil. CONCLUSIONS The spread pattern of DENV-4 Genotype IIB strains across the country supports two independent introductions of the virus in São Paulo in a short period of time. Furthermore, São Paulo appears to have played a pivotal role in the dissemination of DENV-4 to other locations in Brazil

Meta-transcriptomic identification of Trypanosoma spp. in native wildlife species from Australia

Background: Wildlife species carry a remarkable diversity of trypanosomes. The detection of trypanosome infection in native Australian fauna is central to understanding their diversity and host-parasite associations. The implementation of total RNA sequencing (meta-transcriptomics) in trypanosome surveillance and diagnosis provides a powerful methodological approach to better understand the host species distribution of this important group of parasites. Methods: We implemented a meta-transcriptomic approach to detect trypanosomes in a variety of tissues (brain, liver, lung, skin, gonads) sampled from native Australian wildlife, comprising four marsupials (koala, Phascolarctos cinereus; southern brown bandicoot, Isoodon obesulus; swamp wallaby, Wallabia bicolor; bare-nosed wombat, Vombatus ursinus), one bird (regent honeyeater, Anthochaera phrygia) and one amphibian (eastern dwarf tree frog, Litoria fallax). Samples corresponded to both clinically healthy and diseased individuals. Sequencing reads were de novo assembled into contigs and annotated. The evolutionary relationships among the trypanosomatid sequences identified were determined through phylogenetic analysis of 18S rRNA sequences. Results: We detected trypanosome sequences in all six species of vertebrates sampled, with positive samples in multiple organs and tissues confirmed by PCR. Phylogenetic analysis indicated that the trypanosomes infecting marsupials were related to those previously detected in placental and marsupial mammals, while the trypanosome in the regent honeyeater grouped with avian trypanosomes. In contrast, we provide the first evidence for a trypanosome in the eastern dwarf tree frog that was phylogenetically distinct from those described in other amphibians. Conclusions: To our knowledge, this is the first meta-transcriptomic analysis of trypanosomes in native Australian wildlife, expanding the known genetic diversity of these important parasites. We demonstrated that RNA sequencing is sufficiently sensitive to detect low numbers of Trypanosoma transcripts and from diverse hosts and tissues types, thereby representing an effective means to detect trypanosomes that are divergent in genome sequence